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^; ABSTRACT 

O ' 

\^ ■ Aims. 'Water in Star-forming regions with HerscheV (WISH) is a Herscliel Key Programme aimed at understanding the physical and chemical 
"tji structure of young stellar objects (YSOs) with a focus on water and related species, 
j^ Methods. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herscliel Space 

1 I Observatory to measure emission in HtO, CO, OH, [Oi], and [Cn] lines located between 63 and 186 yum. The excitation and spatial distribution 

of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. 

^~^ Results. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the 

^ central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, 

^~^ [Oi] emission is seen in low- velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km s"' in both the red- and 

CO blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical 

^^ envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels 

Cn at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the 

1<~, envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [O i] emission is attributed to /-type shocks in 

^— «) dense gas close to the protostar. In the scenario described here, the combined cooling by far-lR lines within the central spatial pixel is estimated to 

^— «) be 2x10"^ Lq, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope. 

^~~^ Key words. Astrochemistry — Stars: formation — ISM: molecules — ISM: jets and outflows — ISM: individual objects: HH 46 

> 

rS ' 1. Introduction the physical and chemical characteristics ofyoung stellar objects 

rrt (YSOs) as functions of evolutionary stage and across a wide 

. >y The embedded phase of star formation is a critical period m the ^^ luminosities and masses (see also Nisini et al. 2010; 

evolution or a young star, because it is the stage where the hnal p- . j. ^ ooiO) 

mass of the star, the size and mass of the protoplanetary disk, ^i^i • r , t r ^ c^ .^i. 

,,...,, ., .. -..J.,. .. One of the many surprises of the Infrared Space Observatory 

and the initial chemical composition of the disk are determined ^t^^x . j • ^.11 ■. j^^ jtt^i- 

(Andre et al. 2000; Visser et al. 2009; J0rgensen et al. 2009). ^^^^J 7"' *" '^"t'cA"" u ? r^ """^ w , u ? 

-., .1 • 1^ 1 • ^u • J- wards low-mass YSOs with the Long Wavelength Spectrometer 

Many physical processes occur simultaneously in the immedi- y,„rc, ^. ■ ■ , ^r^r^r^ ^ ,,• , i^^r^r^ xt- • • 

. J- f .u . . • f 11 • .u 11 • (LWS; e.g., Giannim et al. 1999; Ceccarelh et al. 1999; Nisini 

ate surroundings of the protostar: inlall in the collapsing enve- , '„„%, „, . . -■ , • , , ^^^ ^^^^ T^^ 1 1 

, .a ^ . ^ J u 1 • .u . • 1 J et al. 2002). The origin of this hot gas (-500-2000 K) has been 

lope, outflows sweeping up and shocking the material, and ener- , ., , , , •■ . ^ ' . „ , ,., 

^- /TT17 J A^ \ u . u .• J J- • .• .u heavily debated, with two main explanations put forward: (1) 

getic (UV and X-ray) photons heating and dissociating the gas , , ,. , , . , f . /-n , 



shocks extending over a large area (arcmin scale), and (ii) the en- 
velope heated by the protostellar luminosity on scales of 2 -20". 
The ISO-LWS data and subsequent modelling could not distin- 



(Spaans et al. 1995; Bachiller & Tafafla 1999; Arce et al. 2007). 

Because of high extinction, these processes can only be probed 

at far-infrared and millimetre wavelengths, but lack of observa- . , , , • • 1 1 1 ^ nn,, 

, c v.- u u J .u • . -c .• T-u ^ c guish bctwccn these scenarios given the large beam of ~80 . 

tional facilities has hampered their quantmcation. The goal of 5,, ,, ^ ., , 7r..^r, n • • .c . 

.u iWT . T o. r • • -.u TT 1 i> /xiiToiix 1 Thc smallcr apcrturc of //er5cne/-PACS allows imaging of thcsc 

the Water In Star-forming regions with //ericne/ (WISH) key ,. r^, a 1^ ■ , , j^ 

■ . ■! r^ /-^A J 1 . J • . J . • lines at 9.4, sufticient to separate the on- and ott-source emis- 

programme is to use HiO, CO and related species to determine . ,, 11 ^ . . . .,-, , , 

^ ° ^ sion. Moreover, its higher sensitivity and higher spectral resolu- 

* Herschel is an ESA space observatory with science instruments "^1°" ^^^°^ detection of weaker lines, 
provided by European-led Principal Investigator consortia and with im- This letter presents Science Demonstration Phase observa- 

portant participation from NASA. tions of HH 46, an isolated low-mass protostar (Lboi ~ 16 Lq, 
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Table 1. Line fluxes observed towards HH 46 with PACS. 
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Fig. 1. The [Oi] 63 /im, CO 16-15 163 yum, H2O 2i2-loi 179 yum lines 
and the OH doublet at 84 fira (possibly blended with a small contribu- 
tion from CO 31-30), at the central position (08'"25™43^:9; -5rOO'36"; 
J2000) and integrated over the red- and blue-shifted outflow lobes. The 
vertical lines show the velocity of [O i] emission at rest (dashed) and 
in the blue- and red-shifted jets (dotted). The Spitzer image of HH 46 
(Velusamy et al. 2007) is shown with the PACS footprint and ISO beam 
overlaid. The spaxels covering the red- and blue-shifted lobes are indi- 
cated with coloured boxes. The central spaxel is shown in black. 



D X 450 pc) located in a dense core with a prominent out- 
flow extending out to ~2'. It has been imaged with the Spitzer 
Space Telescope (Velusamy et al. 2007) and ground-based sub- 
millimetre telescopes (van Kempen et al. 2009, and references 
therein). Because of its well-defined geometry, HH 46 provides 
an ideal testbed for separating the different physical components 
contributing to the observed emission and for benchmarking new 
models of YSOs that can be used for other sources. 



2. Observations and results 

HH 46 was observed on 26 October 2009 with Herschel (Pilbratt 
et al. 2010) using PACS (Poglitsch et al. 2010) in pointed Une- 
scan spectroscopy mode (obsid 1342186315 and 1342186316). 
PACS is a 5x5 aiTay of 9'.'4x9'.'4 spatial pixels (spaxels) that 
cover the 53-210 yum wavelength range with a spectral resolu- 
tion ranging from 1000 to 4000 (the latter only at ~63 //m) in 
spectroscopy mode. In one exposure, a wavelength segment is 
observed in the first order (105-210//m) and at the same time in 
the second (72-105 yum) or third order (53-72 yum). The PACS 
spectrum of HH 46 covers 27 segments obtained in 15 separate 
integrations of 350 s each. Two different nod positions, located 
6' in opposite directions from the target, were used to correct for 
telescopic background. 

Data were reduced with HIPE v2.4.0. The relative spectral 
response function within each band was determined from ground 
calibration prior to launch. The absolute wavelength scale is 
accurate to 30-50 km s"', depending on the position of the 
emission peak in the cross-dispersion axis within the slit. The 
absolute flux calibration below and above 100 yum was sepa- 
rately determined from in-flight observations of (point) calibra- 
tion sources. The uncertainty in absolute and relative fluxes is 
estimated to be 50%, based on a comparison with the ISO-LWS 
continuum data from Nisini et al. (2002). 
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Blue-shifted [O i] component at 



'^ W m"^. Values with colons indicate weak detections 
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Emission is detected in lines of H2O, CO, OH, O i, and C 11 
(Table 1, Fig. 1). Nisini et al. (2002) detected only the [Oi] 63 
and 145 /im and [Cii] 158 /im lines with ISO. The typical 3cr 
sensitivity is ~10"'^ W m"^ to an unresolved emission line and 
unlike the [O i] line, the molecular lines are all unresolved. The 
PACS emission in most lines and in the continuum is strongly 
peaked at the central position. Emission in many lines is also 
seen along the red and blue outflow lobes including bright emis- 
sion in the spaxel centred 11" SW of the source. Tables 1 and 
2 (the latter available online) list the line fluxes at the source 
position, in the red and blue outflows, and in the total field-of- 
view. The fluxes in the central spaxel were corrected for the 
point-source PSF. Fluxes for the outflows were measured over 
the spaxels indicated in Fig. 1 and corrected for the leaking of 
light from the central spaxel into adjacent spaxels. 

Nine lines of water are detected in the central spaxel with 
£'u/A;b ranging from 1 14 to 320 K. Most H2O emission lines peak 
at the location of strong outflow emission 11" SW (Table 2). 
Emission in CO lines ranging from J^ - 14 to 30 (Eu/ks = 580- 
2600 K) is seen in both the central spaxel and in the outflows 
(Tables 1 and 2; Fig. 2). Strong OH emission is detected in the 
four doublets at 79, 84, 119, and 163 /zm, arising from levels up 
to 290 K above the ground state. The OH emission is strongly 
concentrated on-source, although emission in the 119 fim dou- 
blet is also seen in the direction of the red outflow. 

Emission in the [O i] 63 yum and [C 11] 158 fim lines is found 
to be extended over most of the PACS field-of-view. The total 
flux in these two lines is about 1.5 and 7 times weaker, respec- 
tively, in the PACS field (50"x50") than it was in the ISO-LWS 
beam (~80"). Thus, most of the [Cii] and some of the [Oi] 
emission must be located within the ISO field but beyond that 
of PACS or in a relatively smooth background extending out to 
the PACS nod positions. 

The high PACS spectral resolution at 63 fim of ~ 100 km s"' 
allows the [O i] 63 fim emission to be resolved into three veloc- 
ity components (see Fig. 1). In addition to the main component 
around u ^ Q km s"' peaking on-source, strong red- and blue- 
shifted [Oi] 63 yum emission is detected at 110 and -170 km s"', 
respectively, in several spaxels to the SW and NE of the central 
source. Velocities are consistent with the jet velocities measured 
in near-IR and optical lines (e.g., Nishikawa et al. 2008; Garcia 
Lopez et al. 2010). 



3. Analysis 
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Previous observations of HH 46 have revealed many properties 
of the circumstellar environment including the dense centrally 
concentrated envelope containing the protostar HH 46 IRS, the 
presence of warm gas along the outflow walls heated by UV ra- 
diation (van Kempen et al. 2009), the shape and size of the out- 
flow cavities, and the presence of jets and shocks (e.g., Velusamy 
et al. 2007). The temperature (10-250 K) and density (-lO'^-lO'^ 
cm"-*) structure of the passively heated envelope have been de- 
termined in spherical symmetry by fitting dust radiative transfer 
models to the spectral enenergy distribution of the source and the 
spatial extent of the continuum emission. Using these properties, 
a set of existing models (Kristensen et al. 2008; van Kempen 
et al. 2009; Bruderer et al. 2009) is adapted to predict the emis- 
sion in lines of CO, H2O, OH, and O i. In the following, focus is 
placed on analysing emission from the central spaxel. 

Line emission is expected to originate in the known circum- 
stellar components: the passively heated spherical envelope, the 
UV-heated cavity walls, the small-scale shocks along the cavity 
walls, the jet, and the disk. The jet-driven 7-type shock is dis- 
cussed in Sects. 3.3 and 3.4. Based on the PACS observations of 
the HD 100546 disk (Sturm et al. 2010), any disk contribution is 
likely negligible at the distance of HH 46; hence, only the first 
three components are expected to cause the molecular emission. 

For the spherical envelope, the model of van Kempen et al. 
(2009) is rerun with the new 3D non-LTE radiative transfer code 
LIME (Brinch & Hogerheijde in prep.) to obtain fluxes of the 
higher-7 CO lines. The second component, the UV-heated gas in 
the outflow cavity walls (Spaans et al. 1995; van Kempen et al. 
2009), is modelled according to the method of Bruderer et al. 
(2009). The basis is the same spherical envelope profile, but a 
65 000 AU X 13 000 AU ellipsoidal cavity is now carved out 
at an inclination of 53° (Velusamy et al. 2007; Nishikawa et al. 
2008). The only free parameter in this scenario is the protostellar 
FUV luminosity, which is assumed to be 0.1 L© (i.e.. Go ~ lO'* 
at 100 AU with respect to the interstellar radiation field). The 
gas temperature in the cavity walls is parameterised from the 
grid of PDR models by Kaufman et al. (1999) and is typically a 
few hundred K; the dust temperature and density profiles are un- 
changed from the spherical model. More details will be reported 
by Visser et al. (in prep.), who will explore a wider parameter 
space to assess the viability of other scenarios. 

Small-scale shocks created by the outflow along the cav- 
ity walls are the third component considered for the molecular 
emission. Their temperature is typically a few thousand K. The 
shock emission is computed by tiling a number of ID C-type 
shock models along the 2D cavity shape (Kristensen et al. 2008), 
taking the width of each shock to be the region over which the 
considered species contributes significantly to the cooling in ID 
shock models. This can elTectively be considered as an estimate 
of the beam filling factor. For each density in the range lO'^-lO^^ 
cm"^, the emission is computed using the results from Kaufman 
& Neufeld (1996). The only free parameter in this model is the 
shock velocity, which is assumed constant along the walls. For 
the case of CO, a velocity of 20 km s ' reproduces the observa- 
tions, and this velocity is adopted for the other species as well. 
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Fig. 2. CO line fluxes observed in the central PACS spaxel (J^ > 10) 
and with APEX (/„ < 10). Model fluxes are used to estimate the ratio 
of flux in a Active PACS spaxel at the APEX wavelength and the ob- 
served APEX flux. Overplotted are predictions from a passively heated 
envelope (blue), a UV-heated cavity (green), and small-scale shocks in 
the cavity walls (red). The black line is the sum of the three. A cartoon 
of the different components is shown in the inset. 



and below 10^ cm"^. In colder regions, freeze-out lowers the 
CO abundance by a factor of 100. 

The model spectra are convolved with the PACS beam at the 
relevant wavelength and extracted from a 9'.'4 square spaxel at 
the centre. Figure 2 shows the observations, together with the 
model predictions from the passive envelope, the UV-heated cav- 
ity walls, and the small-scale C-type shocks. Individually, each 
component only fits part of the data, but together they reproduce 
the observations over the entire range of rotational levels from 
Ju - 2-32. The results confirm the plausibility of the scenario 
without excluding other solutions not investigated here. 

3.2. H2O 

The passive envelope underproduces the observed H2O fluxes by 
two orders of magnitude. Predicting fluxes from any 2D model 
such as the UV-heated cavity walls is uncertain by an order 
of magnitude due to challenges of radiative transfer modelling 
of H2O. Within these uncertainties, both the UV-heated cavity 
model and the C-type shock model are able to reproduce the ob- 
servations independently. The former requires an H2O gas abun- 
dance of only ~10"^ in the cavity walls and ~10"** in the rest 
of the envelope, as expected from chemical models including 
photodissociation and freeze-out. The C-type shock component 
matches the observations if the abundances from Kaufman & 
Neufeld (1996) are scaled down to ~ 7x10"^. This could be 
accomplished by photodissociation of H2O in the shocked gas. 
More detailed modelling, including spatially extended emission 
and a comparison with spectrally resolved line profiles observed 
with HIFI, is needed to distinguish these scenarios. 



3.1. CO 

The PACS data (Sect. 2) are complemented by spectrally re- 
solved /u < 7 lines (E^/kB < 155 K) detected with APEX 
(Av a 10 km s '; van Kempen et al. 2009). The CO/H2 abun- 
dance ratio in the model is taken to be 1.6x10"'* above 20 K 



3.3. [O I] 

Within the passively heated envelope and cavity walls, the 2D 
models presented above yield an [O i] 63 fim emission line that 
is narrow (3^ km s"'), optically thick, and weaker than the ob- 
served flux by more than three orders of magnitude, for an O i 
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abundance of 10 *". C-type shocks fail to reproduce the obser- 
vations by seven orders of magnitude, and are still insufficient 
even if all H2O is dissociated into O i. The observed emission 
is, however, reproduced in a dissociative /-type shock (Neufeld 
& Dalgarno 1989). Thus, for the central spaxel, where bright 
emission is seen at the systemic velocity, [Oi] may trace the 
impact of the high-velocity jet on the densest parts of the en- 
velope (riu > 10^ cm"-'), causing rapid deceleration from u > 
200 km s ' while cooling the gas through [Oi] emission. The 
observed [Oi] 63/145 fim line ratio is ~16 in the central spaxel. 
Shock models predict ratios of 14-20, where the lower ratio ap- 
plies to higher densities (10*" cm"^; Neufeld & Dalgarno 1989). 
The high-y [Oi] emission observed in both the central and adja- 
cent spaxels can be produced in fast, dissociative shocks in the 
much lower density jet itself. 



Table 3. Origin of line emission and cooling rates of various species for 
the HH 46 central spaxel. 



Species 



Total 



Cooling 

10-3 u, 



Origin 



CO 


0.1 


CO 


3.8 


CO 


2.8 


H20 


0.1 


H20 


5.0 


Oi 


9.5 


OH 


>2.4 


Cii 


0.1 



Passively heated envelope 
UV-heated cavity walls 
C-type shocks 
Passively heated envelope 
UV or C-type shocks 
y-type shock 
y-type shock 
Surrounding cloud 



>23.8 



3.4. OH 

The OH line ratios were modelled separately for a single- 
component slab model using an escape probability code with 
absorbing and emitting dust continuum (similar to Bruderer et 
al. subm.). Comparison with the observed line ratios and inten- 
sities shows that the OH lines likely originate in a high-density 
and high-te mperature region (nH~ 10^ cm"^, rgas>800K). Along 
with the lack of extended emission, this rules out an origin in the 
photon-heated cavity walls and C-type shocks along the cavity 
walls. An OH column density of ~10'^ cm"^ gives the best fit to 
the observations along with a physical size of the OH emitting 
region of O'.'S (~250 AU). Other solutions are possible and will 
be discussed in a forthcoming paper (Wampfler et al. in prep.). 

The models indicate that better fits are obtained for tempera- 
tures higher than what can be accounted for in a passive envelope 
model. The major competing coolant to Oi in a fast, dissocia- 
tive shock is OH (Neufeld & Dalgarno 1989), and it is therefore 
likely that some of the OH emission is caused by the jet shock 
impinging on the inner, dense envelope. Further modelling, in- 
cluding combinations of C- and /-type shocks (Snell et al. 2005), 
will be explored to constrain the origin of the OH emission. 



3.5. Origin of hot gas 

Table 3 summarizes the assigment of the various species to the 
different proposed physical components. It also includes the total 
far-IR cooling through H2O, CO, OH, and [O i] lines in the cen- 
tral spaxel. Within our scenario, cooling by CO takes place al- 
most equally through the UV-heated cavity walls and the small- 
scale C-type shocks. The observed H2O emission can be ac- 
counted for either by UV heating of the cavity walls or by C- 
type shocks, or a combination. It is impossible to distinguish 
between the two scenarios at present; however, the total cooling 
is SxlO""' Lq in both cases using the above models to account 
for the non-observed line emission. 

Cooling through [Oi] emission takes place on very small 
spatial scales, probably related directly to the jet impinging on 
the envelope walls. OH emission likely arises in the same disso- 
ciative shock. The OH cooling in observed lines is 2.4x10"^ L©, 
but the total OH cooling can be an order of magnitude higher if 
an excitation temperature of 200 K is used to account for non- 
observed lines. Thus, the total cooling caused by J- and C-type 
shocks is at least 1.5x10"^ Lq or 60% of the total far-lR line 
cooling. The two other components (passive and UV) are re- 
sponsible for the remaining 40% of the total cooling. ISO-LWS 
found typical far-lR line cooling rates of 1-5x10"^ L© (Nisini 



et al. 2002), comparable to what is observed here. The total cool- 
ing of the entire HH 46 system is higher, since only the central 
spaxel is considered here. 

In summary, the Herschel-PACS data allow for disentangling 
and quantifying the energetic processes occurring in deeply em- 
bedded protostars. The models indicate, that of the two scenarios 
previously proposed based on ISO-LWS data, shocks are more 
important than the passively heated envelope in powering the far- 
infrared lines. Another component, the UV-heated cavity walls, 
also plays an important role in producing line emission. HlFl 
observations resolving emission lines are planned to test the pro- 
posed scenario in the framework of the WISH key programme. 

Acknowledgements. This work is made possible thanks to the HIFI guaranteed 
time programme and the PACS instrument builders. We thank Javier Gracia 
Carpio, Jeroen Bouwman, Bruno Merin, and Bart VandenBussche for help with 
the PACS data reduction and many funding agencies for financial support. We 
also thank the referee for comments that improved this paper 



References 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, in Protostars & Planets IV, 

ed. V. Mannings, A. P Boss, & S. S. Russell (Tucson: UAP), 59 
Arce, H. G., Shepherd, D., Gueth, R, et al. 2007, in Protostars & Planets V, ed. 

B. Reipurth, D. Jewitt, & K. Keil (Tucson: UAP), 245 
Bachiller, R. & Tafalla, M. 1999, in The Origin of Stars and Planetary Systems, 

ed. C. J. Lada, , & N. D. Kylafis (Kluwer, Dordrecht), 227 
Bruderer, S., Benz, A. O., Doty, S. D., van Dishoeck, E. F., & Bourke, T. L. 

2009, ApJ, 700, 872 
CeccareUi, C, Caux, E., Loinard, L., et al. 1999, A&A, 342, L21 
Fich et al. 2010, A&A, this volume 

Garcia Lopez, R., Nisini, B., Eisloflfel, J., et al. 2010, A&A, 511, A5 
Giannini, T, Lorenzetti, D., Tommasi, E., et al. 1999, A&A, 346, 617 
J0rgensen, J. K., van Dishoeck, E. P., Visser, R., et al. 2009, A&A, 507, 861 
Kaufman, M. J. & Neufeld, D. A. 1996, ApJ, 456, 611 
Kaufman, M. J., Wolfire, M. G., HoUenbach, D. J., & Luhman, M. L. 1999, ApJ, 

527, 795 
Kristensen, L. E., Ravkilde, T. L., Pineau des Forets, G., et al. 2008, A&A, 477, 

203 
Neufeld, D. A. & Dalgarno, A. 1989, ApJ, 344, 251 
Nishikawa, T, Takami, M., Hayashi, M., Wiseman, J., & Pyo, T 2008, ApJ, 684, 

1260 
Nisini, B., Giannini, T, & Lorenzetti, D. 2002, ApJ, 574, 246 
Nisini et al. 2010, A&A, this volume 
Pilbratt et al. 2010, A&A, this volume 
Poglitsch et al. 2010, A&A, this volume 

Snell, R. L., HoUenbach, D., Howe, J. E., et al. 2005, ApJ, 620, 758 
Spaans, M., Hogerheijde, M. R., Mundy, L. G., & van Dishoeck, E. F. 1995, 

ApJ, 455, 167 
Sturm et al. 2010, A&A, this volume 

van Kempen, T A., van Dishoeck, E. F, Glisten, R., et al. 2009, A&A, 501, 633 
Velusamy, T, Langer, W. D., & Marsh, K. A. 2007, ApJ, 668, L159 
Visser, R., van Dishoeck, E. R, Doty, S. D., & Dullemond, C. P 2009, A&A, 

495, 881 



T.A. van Kempen et al.: Herschel-PACS spectroscopy of HH 46 



' Leiden Observatory, Leiden University, PO Box 9513, 2300 RA 
Leiden, Tiie Netherlands 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 
MS 42, Cambridge, MA 02138, USA 

' Max-Planck-Institut fiir Extraterrestrische Physik, 

Giessenbachstrasse 1, 85748 Garching, Germany 

"* Institute of Astronomy, ETH Zurich, 8093 Zurich, Switzerland 

^ Department of Physics and Astronomy, Denison University, 
Granville, OH, 43023, USA 

* Centre for Star and Planet Formation, Natural History Museum of 
Denmark, University of Copenhagen, 0ster Voldgade 5-7, DK-1350 
Copenhagen K., Denmark 

' Observatorio Astronomico Nacional (IGN), Calle Alfonso Xn,3. 
28014, Madrid, Spain 

** Department of Physics and Astronomy, Johns Hopkins University, 
3400 North Charles Street, Baltimore, MD 21218, USA 

' Universite de Bordeaux, Laboratoire d'Astrophysique de 
Bordeaux, France; CNRS/INSU, UMR 5804, Floirac, France 
"^ INAF - Instituto di Fisica dello Spazio Interplanetario, Area di 
Ricerca di Tor Vergata, via Fosso del Cavaliere 100, 00133 Roma, 
Italy 

" Department of Astronomy, University of Michigan, 500 Church 
Street, Ann Arbor, MI 48109-1042, USA 

'^ Department of Radio and Space Science, Chalmers University of 
Technology, Onsala Space Observatory, 439 92 Onsala, Sweden 
'^ California Institute of Technology, Division of Geological and 
Planetary Sciences, MS 150-21, Pasadena, CA 91125, USA 
''' School of Physics and Astronomy, University of Leeds, Leeds 
LS2 9JT, UK 

'^ INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 
50125 Firenze, Italy 

'* Centro de Astrobiologia, Departamento de Astrofisica, CSIC- 
INTA, Carretera de Ajalvir, Km 4, Torrejon de Ardoz. 28850, 
Madrid, Spain 

" Astronomical Institute Anton Pannekoek, University of 
Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 
"* Department of Astrophysics/IMAPP, Radboud University 
Nijmegen, PO. Box 9010, 6500 GL Nijmegen, The Netherlands 
" LERMA and UMR 81 12 du CNRS, Observatoire de Paris, 61 Av. 
de r Observatoire, 75014 Paris, France 

^"^ University of Waterloo, Department of Physics and Astronomy, 
Waterloo, Ontario, Canada 

^' Observatorio Astronomico Nacional, Apartado 1 12, 28803 Alcala 
de Henares, Spain 

^^ INAF - Osservatorio Astronomico di Roma, 00040 Monte Porzio 
catone, Italy 

^' SRON Netherlands Institute for Space Research, PO Box 800, 
9700 AV, Groningen, The Netherlands 

^"' National Research Council Canada, Herzberg Institute of 
Astrophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, 
Canada 

^' Department of Physics and Astronomy, University of Victoria, 
Victoria, BC V8P lAI, Canada 

^^ Department of Physics and Astronomy, San Jose State University, 
One Washington Square, San Jose, CA 95192, USA 
^^ Department of Astronomy, Stockholm University, AlbaNova, 106 
91 Stockholm, Sweden 

^^ California Institute of Technology, Cahill Center for Astronomy 
and Astrophysics, MS 301-17, Pasadena, CA 91125, USA 
^' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 
53121 Bonn, Germany 

^^ Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, CA 91109, USA 

^' Department of Physics and Astronomy, University of Calgary, 
Calgary, T2N 1N4, AB, Canada 

^^ Instituto de RadioAstronomia Milimetrica, Avenida Divina 
Pastora, 7, Nucleo Central E 18012 Granada, Spain 
^^ Kapteyn Astronomical Institute, University of Groningen, PO 
Box 800, 9700 AV, Groningen, The Netherlands 



'** KOSMA, I. Physik. Institut, Universitat zu Koln, Zulpicher Str. 

77, D 50937 Koln, Germany 

35 California Institute of Technology, 1200 E. California Bl., MC 

100-22, Pasadena, CA. 91125 USA 

3* Experimental Physics Dept., National University of Ireland 

Maynooth, Co. Kildare. Ireland 



T.A. van Kempen et al.: Herschel-PACS spectroscopy of HH 46, Online Material p 1 

Table 2. Line fluxes observed towards HH 46 with PACS." 

Central Spaxel Red Outflow Blue Outflow All 25 Spaxels 

Area (square arcsec) 88 530 530 2200 

Species Transition /Ij^b E^/kg Flux Flux Flux Flux 

Qjm) (K) (IQ-'^Wm-^) 

~CO 14^43 185.999 5805 48^4 27T2 r8T2 108 ±7 

16-15 162.812 751.7 46 ± 5 27 ± 2 13 + 3 80 ± 8'' 

18-17 144.784 945.0 71 ± 7 27+2 <14 67 ± 30'' 

22-21 118.581 1397.4 33 ± 7 <34 <27 44 ± 8 

23-22*= 113.458 1524.2 78 ±11 15 ± 3 <14 78 ± 13'' 





24-23 


108.763 


1656.5 


18 ±6 


<14 


<26 


38 ±9 




29-28 


90.163 


2399.8 


18±3 


<26 


<19 


28 ±6 




30-29 


87.190 


2564.8 


25 ±4: 


<14 


<21 


<55 




31-30'' 


84.411 


2735.3 


74 ±6 


<31 


<22 


55 ±7" 




32-31 


81.806 


2911.2 


<29 


<40 


<40 


<61 




33-32 


79.360 


3092.5 


<17 


<27 


<41 


<57 




36-35 


72.843 


3668.8 


<23 


<49 


<51 


<136 


H20 


221-2,2 


180.488 


194.1 


<7 


<9 


<9 


<14 




2i2-loi 


179.527 


114.4 


4 ±2: 


46 ±3 


16 ±2 


82 ±7 




3o3-2l2 


174.626 


196.8 


<22 


24 ±2 


<21 


54 ±7 




3l3-2o2 


138.528 


204.7 


14 ±4 


18 ±3: 


<10 


34 ±8 




4o4-3l3 


125.354 


319.5 


<18 


<6 


14 ±3 


<24 




4i4-3o3'-' 


113.537 


323.5 


78 ±11 


15 ±3 


<14 


78 ± 13" 




221-1 10 


108.073 


194.1 


17 ±7 


16 ±3: 


21 ±4 


64 ±11 




322-211 


89.988 


296.8 


<27 


<22 


<49 


<110 




8l8-7o7 


63.324 


1070.7 


<33 


<10 


<22 


<43 


Off 


13 11 


/163.396 
163.015 


269.8 


<27 


<4 


<9 


<10 




2'2 2'2 


270.2 


22 ±4 


<18 


<5 


20 ±4 




3 5 3 3 


f 119.441 
119.234 


120.5 


44 ±7 


49 ±5 


<7 


86 ± 12'' 




2'2 2'2 


120.7 


38 ±9 


25 ±4 


<19 


71 ±15 




3 7 3 5d 


f 84.597 
84.420 


290.5 


87 ±6 


<14 


<27 


95 ±10 




2'2 2'2 


291.2 


74 ±6 


<31 


<22 


55 ±7" 




113 3 


/ 79.179 
79.116 


181.7 


38 ±5 


<41 


<37 


93 ±13 




2'2 2'2 


181.9 


55 ±7 


<34 


<35 


87 ±10 


Oi 


'Po-^P2 


145.525 


326.6 


82 ±8 


68 ±4 


48 ±3 


208 ± 14 




3P1-3P2 


63.184 


227.7 


1260 ±54 


186 ±28 


420 ± 20 


1870 ± 80 






!) = -170kms-' f 


162 ±12 


- 


90 ±8 


282 ± 51 






!; = +110kms-'f 


- 


1230 ± 52 


- 


1510 ±100 


Cii 


'P3/2-'Pl/2 


157.741 


91.2 


73 ±6 


77 ±4 


140 ±6 


505 ± 23 



Fluxes are measured from Gaussian fits. Lines marked with a colon are only weakly detected. Listed uncertainties are 68% confidence intervals, 

upper limits are 95% confidence intervals. The uncertainties do not include the 50% error margin in the relative spectral response function. 

The total flux is less than the combined flux from the central spaxel, the red outflow and the blue outflow. This is a result of higher noise in the 

spectrum summed over all spaxels, leading to a different fit and PSF over-correction. 

The CO 23-22 and H2O 4i4-3o3 lines at 113.5 //m are blended. 

The CO 31-30 and OH 3/2,7/2-3/2,5/2 lines at 84.4 yum are blended. 

Doublet transitions. 

Blue- and red-shifted components of the [O i] 63.2 /jm line. The blue component is not detected in the spaxels covering the red outflow, and 

vice versa. Futhermore, the red component is not detected in the central spaxel. 



